(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     27354,        722]
NotebookOptionsPosition[     26100,        674]
NotebookOutlinePosition[     26437,        689]
CellTagsIndexPosition[     26394,        686]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{"Integrate", "[", 
  RowBox[{
   RowBox[{"Log", "[", 
    RowBox[{"Sqrt", "[", 
     RowBox[{"2", "-", 
      RowBox[{"2", " ", 
       RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], "]"}], "]"}], ",", 
   "\[Theta]"}], "]"}]], "Input",
 CellChangeTimes->{{3.515419091485246*^9, 3.515419092074631*^9}, {
   3.515419152981728*^9, 3.515419153750572*^9}, {3.515419377245436*^9, 
   3.515419458950328*^9}, {3.515419547941479*^9, 3.51541958216004*^9}, {
   3.515419672726908*^9, 3.515419675631078*^9}, {3.515419736852989*^9, 
   3.515419743787541*^9}, {3.515419820801923*^9, 3.515419853242006*^9}, {
   3.515419934279351*^9, 3.515419971575825*^9}, 3.515420020212903*^9}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "2"], " ", 
  RowBox[{"(", 
   RowBox[{
    FractionBox[
     RowBox[{"\[ImaginaryI]", " ", 
      SuperscriptBox["\[Theta]", "2"]}], "2"], "-", 
    RowBox[{"2", " ", "\[Theta]", " ", 
     RowBox[{"Log", "[", 
      RowBox[{"1", "-", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"\[ImaginaryI]", " ", "\[Theta]"}]]}], "]"}]}], "+", 
    RowBox[{"\[Theta]", " ", 
     RowBox[{"Log", "[", 
      RowBox[{"2", "-", 
       RowBox[{"2", " ", 
        RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], "]"}]}], "+", 
    RowBox[{"2", " ", "\[ImaginaryI]", " ", 
     RowBox[{"PolyLog", "[", 
      RowBox[{"2", ",", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"\[ImaginaryI]", " ", "\[Theta]"}]]}], "]"}]}]}], 
   ")"}]}]], "Output",
 CellChangeTimes->{3.515420020730261*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"F", "[", "\[Theta]_", "]"}], ":=", " ", 
   RowBox[{
    FractionBox["1", "2"], " ", 
    RowBox[{"(", 
     RowBox[{
      FractionBox[
       RowBox[{"\[ImaginaryI]", " ", 
        SuperscriptBox["\[Theta]", "2"]}], "2"], "-", 
      RowBox[{"2", " ", "\[Theta]", " ", 
       RowBox[{"Log", "[", 
        RowBox[{"1", "-", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Theta]"}]]}], "]"}]}], "+", 
      RowBox[{"\[Theta]", " ", 
       RowBox[{"Log", "[", 
        RowBox[{"2", "-", 
         RowBox[{"2", " ", 
          RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], "]"}]}], "+", 
      RowBox[{"2", " ", "\[ImaginaryI]", " ", 
       RowBox[{"PolyLog", "[", 
        RowBox[{"2", ",", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Theta]"}]]}], "]"}]}]}], ")"}]}]}],
   ";"}]], "Input",
 CellChangeTimes->{{3.515419976087289*^9, 3.515419976823832*^9}, {
  3.515420023851779*^9, 3.515420036455968*^9}, {3.515497647666007*^9, 
  3.515497648652904*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"k", "=", 
  RowBox[{"FullSimplify", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"n", "^", "2"}], "/", 
     RowBox[{"(", "Pi", ")"}]}], "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"F", "[", "Pi", "]"}], "-", 
      RowBox[{"F", "[", 
       RowBox[{"Pi", "/", "n"}], "]"}]}], ")"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.515419979483451*^9, 3.515420001595596*^9}, {
  3.515420043089247*^9, 3.515420046975592*^9}, {3.515420080520872*^9, 
  3.515420174624087*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   FractionBox["1", "12"], " ", "\[ImaginaryI]", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "3"}], "+", 
     RowBox[{"2", " ", 
      SuperscriptBox["n", "2"]}]}], ")"}], " ", "\[Pi]"}], "+", 
  RowBox[{"n", " ", 
   RowBox[{"Log", "[", 
    RowBox[{"1", "-", 
     SuperscriptBox["\[ExponentialE]", 
      FractionBox[
       RowBox[{"\[ImaginaryI]", " ", "\[Pi]"}], "n"]]}], "]"}]}], "-", 
  RowBox[{
   FractionBox["1", "2"], " ", "n", " ", 
   RowBox[{"Log", "[", 
    RowBox[{"2", "-", 
     RowBox[{"2", " ", 
      RowBox[{"Cos", "[", 
       FractionBox["\[Pi]", "n"], "]"}]}]}], "]"}]}], "-", 
  FractionBox[
   RowBox[{"\[ImaginaryI]", " ", 
    SuperscriptBox["n", "2"], " ", 
    RowBox[{"PolyLog", "[", 
     RowBox[{"2", ",", 
      SuperscriptBox["\[ExponentialE]", 
       FractionBox[
        RowBox[{"\[ImaginaryI]", " ", "\[Pi]"}], "n"]]}], "]"}]}], 
   "\[Pi]"]}]], "Output",
 CellChangeTimes->{{3.515419983311882*^9, 3.515420002076315*^9}, {
  3.515420038498804*^9, 3.515420049866889*^9}, {3.515420087990827*^9, 
  3.515420137153031*^9}, {3.515420168955858*^9, 3.515420177672644*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"k", "-", 
      RowBox[{"n", "*", 
       RowBox[{"(", 
        RowBox[{"1", "-", 
         RowBox[{"Log", "[", "Pi", "]"}]}], ")"}]}]}], ",", 
     RowBox[{"n", "*", 
      RowBox[{"Log", "[", "n", "]"}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"n", ",", "1", ",", "100"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.515420109518128*^9, 3.515420128959294*^9}, 
   3.515420206732399*^9, {3.515420333772683*^9, 3.515420342337107*^9}}],

Cell[BoxData[
 FormBox[
  LineWrapBox["\[AutoOperand]"], InputForm]], "Output",
 CellChangeTimes->{{3.515420120094639*^9, 3.515420139677082*^9}, {
   3.515420170706857*^9, 3.515420207193481*^9}, 3.515420342982319*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Series", "[", 
  RowBox[{"k", ",", 
   RowBox[{"{", 
    RowBox[{"n", ",", "Infinity", ",", "3"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.51542021657657*^9, 3.515420223172666*^9}, {
  3.515420259696464*^9, 3.515420259887567*^9}, {3.51542031704796*^9, 
  3.515420325504303*^9}, {3.515420401094032*^9, 3.515420402155411*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"1", "-", 
      RowBox[{"Log", "[", 
       FractionBox["1", "n"], "]"}], "-", 
      RowBox[{"Log", "[", "\[Pi]", "]"}]}], ")"}], " ", "n"}], "+", 
   FractionBox[
    SuperscriptBox["\[Pi]", "2"], 
    RowBox[{"72", " ", "n"}]], "+", 
   FractionBox[
    SuperscriptBox["\[Pi]", "4"], 
    RowBox[{"14400", " ", 
     SuperscriptBox["n", "3"]}]], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", 
      FractionBox["1", "n"], "]"}], "4"],
    SeriesData[$CellContext`n, 
     DirectedInfinity[1], {}, -1, 4, 1],
    Editable->False]}],
  SeriesData[$CellContext`n, 
   DirectedInfinity[1], {
   1 - Log[$CellContext`n^(-1)] - Log[Pi], 0, Rational[1, 72] Pi^2, 0, 
    Rational[1, 14400] Pi^4}, -1, 4, 1],
  Editable->False]], "Output",
 CellChangeTimes->{3.51542022381215*^9, 3.515420260739414*^9, 
  3.51542032602966*^9, 3.515420403034172*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"H", " ", "=", " ", 
   RowBox[{"FullSimplify", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"n", "^", "2"}], "/", "Pi"}], " ", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"F", "[", "Pi", "]"}], "-", 
       RowBox[{"F", "[", "\[Theta]0", "]"}]}], ")"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.515497620754384*^9, 3.515497642800147*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", 
  RowBox[{
   RowBox[{"Series", "[", 
    RowBox[{"H", ",", 
     RowBox[{"{", 
      RowBox[{"\[Theta]0", ",", "0", ",", "3"}], "}"}]}], "]"}], ",", 
   RowBox[{"Assumptions", "\[Rule]", 
    RowBox[{"\[Theta]0", " ", ">", " ", "0"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.515497660989117*^9, 3.515497697670912*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   RowBox[{"-", 
    FractionBox[
     RowBox[{
      SuperscriptBox["n", "2"], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        RowBox[{"Log", "[", "\[Theta]0", "]"}]}], ")"}], " ", "\[Theta]0"}], 
     "\[Pi]"]}], "+", 
   FractionBox[
    RowBox[{
     SuperscriptBox["n", "2"], " ", 
     SuperscriptBox["\[Theta]0", "3"]}], 
    RowBox[{"72", " ", "\[Pi]"}]], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", "\[Theta]0", "]"}], "4"],
    SeriesData[$CellContext`\[Theta]0, 0, {}, 1, 4, 1],
    Editable->False]}],
  SeriesData[$CellContext`\[Theta]0, 
   0, {-$CellContext`n^2 Pi^(-1) (-1 + Log[$CellContext`\[Theta]0]), 0, 
    Rational[1, 72] $CellContext`n^2/Pi}, 1, 4, 1],
  Editable->False]], "Output",
 CellChangeTimes->{{3.51549767265157*^9, 3.515497700383541*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Reduce", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"-", 
     FractionBox[
      RowBox[{
       SuperscriptBox["n", "2"], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", 
         RowBox[{"Log", "[", "\[Theta]0", "]"}]}], ")"}], " ", "\[Theta]0"}], 
      "\[Pi]"]}], "\[Equal]", 
    RowBox[{"n", "*", 
     RowBox[{"Log", "[", "n", "]"}]}]}], ",", "\[Theta]0", ",", "Reals"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.515497713936037*^9, 3.515497730658506*^9}, {
  3.515497801215959*^9, 3.51549780279486*^9}, {3.515497972399259*^9, 
  3.515497972994489*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"n", "\[Equal]", "1"}], "&&", 
    RowBox[{"\[Theta]0", "\[Equal]", "\[ExponentialE]"}]}], ")"}], "||", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"1", "<", "n", "\[LessEqual]", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"-", 
         RowBox[{"ProductLog", "[", 
          RowBox[{"-", 
           FractionBox["1", "\[Pi]"]}], "]"}]}]]}], "||", 
      RowBox[{"n", "\[GreaterEqual]", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"-", 
         RowBox[{"ProductLog", "[", 
          RowBox[{
           RowBox[{"-", "1"}], ",", 
           RowBox[{"-", 
            FractionBox["1", "\[Pi]"]}]}], "]"}]}]]}]}], ")"}], "&&", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"\[Theta]0", "\[Equal]", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"1", "+", 
         RowBox[{"ProductLog", "[", 
          RowBox[{"-", 
           FractionBox[
            RowBox[{"\[Pi]", " ", 
             RowBox[{"Log", "[", "n", "]"}]}], 
            RowBox[{"\[ExponentialE]", " ", "n"}]]}], "]"}]}]]}], "||", 
      RowBox[{"\[Theta]0", "\[Equal]", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"1", "+", 
         RowBox[{"ProductLog", "[", 
          RowBox[{
           RowBox[{"-", "1"}], ",", 
           RowBox[{"-", 
            FractionBox[
             RowBox[{"\[Pi]", " ", 
              RowBox[{"Log", "[", "n", "]"}]}], 
             RowBox[{"\[ExponentialE]", " ", "n"}]]}]}], "]"}]}]]}]}], 
     ")"}]}], ")"}], "||", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"0", "<", "n", "<", "1"}], "&&", 
    RowBox[{"\[Theta]0", "\[Equal]", 
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{"1", "+", 
       RowBox[{"ProductLog", "[", 
        RowBox[{"-", 
         FractionBox[
          RowBox[{"\[Pi]", " ", 
           RowBox[{"Log", "[", "n", "]"}]}], 
          RowBox[{"\[ExponentialE]", " ", "n"}]]}], "]"}]}]]}]}], 
   ")"}]}]], "Output",
 CellChangeTimes->{
  3.515497731813395*^9, {3.51549796986657*^9, 3.51549797362087*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"g", "[", "n_", "]"}], " ", ":=", 
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{"1", "+", 
     RowBox[{"ProductLog", "[", 
      RowBox[{"-", 
       FractionBox[
        RowBox[{"\[Pi]", " ", 
         RowBox[{"Log", "[", "n", "]"}]}], 
        RowBox[{"\[ExponentialE]", " ", "n"}]]}], "]"}]}]]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"g", "[", "n", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"n", ",", "1", ",", "100"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.515497762921515*^9, 3.515497794140939*^9}, {
  3.515498195490314*^9, 3.515498204819076*^9}, {3.515756840770513*^9, 
  3.515756849551772*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], 
    LineBox[{{1.0000020204081632`, 2.7182754811714283`}, {1.030365074135403, 
     2.6254735781701255`}, {1.0607281278626428`, 2.53752336590795}, {
     1.1214542353171222`, 2.374447430775303}, {1.151817289044362, 
     2.2985590082269725`}, {1.1821803427716018`, 2.2259846380248676`}, {
     1.2429064502260814`, 2.0895795986207073`}, {1.2732695039533213`, 
     2.025199644523845}, {1.303632557680561, 1.9630265784624932`}, {
     1.3643586651350406`, 1.8443287829610642`}, {1.4426195657733265`, 
     1.7014461616150467`}}], LineBox[CompressedData["
1:eJwV0Hk41IkDBvA5vrK2rJwdq2LSMI6pQZgO3zc2lVLO5ot+KxYhIYS2seuX
koQiHaN+lCJklc1VRJoe9CAdjmpLEVq5Ipswx87vj/d5n89/7/Pq+4W7BjBo
NFqUIv/vRqLO06NBDwFFuLyCP0f+2t6eglE9nJ51b6nwmiMfuholLenXR2Ax
p+jQdxJy5w7N7s4pfeSLf204pykhX2+WGWYqsVDa7bogb7mEnLDpePy9IQvO
iZzVCZYSUm9FgqokmAWTHht9Tx8JGT/cndUzzgLv23zV0jIJaXPsxLVrcytR
Nq7zydpRSj6Kj5zyWWCAi3trGfZuUtI59j+bdZcZ4NaVr7bYIyX3hVh8zLI1
gEHrwjL1MCl53vmdcdJ/DUDaZztoZUjJKV3rP4OUV8GsOWS0vENK3qoYfGCm
yUb7RMEbczcZmaNfbNjKYoO5osPc1ktGpqYdSA8xZ6OyROMyfGVkUMCU9w0X
NsI5TpPG4TKSpc2Y1jvDxu4fNty8kCIjz0cvM9NSNYT/m+rU9/UyUmjhLpqd
ZwSHMB3NKwZyMiR3EU2kYwTao0OGpcZy0nP+X4HWbCNkztvkULFGTlp98LWM
3myErFbj7tINcnI8M7x9NNEIBptMbjq7y0m/yRSlXhoHoZOpXZ5H5eS2soaD
TbMcjNn2dT57JSfT0++JyucbozFfqT/5nZzs2H+n4aquMYgpn/38ATm5l52v
LrQ1Rn5Mq/aJz3Iy5lJy2ZqjxrA8eEV9P5OGvOO7JkQqJiCmQ5PYS2iQeL6N
CF1sCtq4dDTKjobbjJlwdSsuqgiPpcFpNBRLzmbccuBCGBVQfuY0Dde/cu84
CbiYCGwmKzNouPDJ/+vJOC6KW1Wt587R8Nvzp0JGDRfc+W+koTk07MgrPPmF
XA33qvsnJ0tpmNIrefbJYA1SBnmOoic0ZDzwHm7l8jBQnSHq/Y6OBzHaymlW
POxz62TGf0/HmGk7y8mWhw8x680WLaDDUWTn1ebEQ55vnniLGh2MSM7jtlAe
RqxKONnadESvnC54UsyDeuwTIyUWHVRSpt9Ttjlistz7pXw69HY8fvViuQVc
SmpyHwbSMXtv2e71hhYwtZnxWhtERwcn8nneagvwRySSgmA6kpV/bDkIC3RG
Z88mhdIx/jC0Vs3PAvbJ6bR1kXTU8Rfmbs+3gMf1Q24H4unwNhT4i00ssTjs
S6r6WTrOMwbG/rRZi0b+uE56DR1+Smr1HVxrXOrpNs6Q0VGeIxRr/8iHoHlt
mI81AyN+F8e4vPUIYdX+T+0AA0eiHnakrNyIrKJDCYtFDHwuze6tI0kob3Wp
CRIzEANvdZUi4JLQNvaXvxnI3RWh5n5lE6Z1nlhGKTOxPvqN3qkSO3yex7jz
zyomGPXtpz8U2GPUrGvwLcnExhVLNYwqf8JgReLYMU8mahyTuMIHm1FulFWS
GcGE/+XDhRPVDlDKy10VfIyJBOajtW6NW1AxtJXtfJGJpjPn9LoatuJqzqFP
aUVMuAwl5wdVb8PZnM/HPe4ysVPwYuNHsSMOB9KGXB8zETN83eLE/e1YPjqP
dOxmQuNc2wbN1h0wtTpaKBxgYkJl8fvSTifcXdLAOzPJxNz1off723aiqW/z
72oyJuI8NP7Rer0LT/c4r+9UIbBQx/dOeIgzbqxexOdrEVD6uO3ungFndCcv
betfRqDcrqe5PsQFtcmFTa8NCcTxrE5YfHQBvTnj8a01BL4W83ctC3OFWBbb
0cIncLKK1rhoxBUn+l1bBjYRqCt6qZK9zw0cFrqE2wiUsPcIMkfdUJf73MfL
hcBtq2tul0LcoZ7Vv86HIhD27Jzn0kl3mAgFJQU+BNpcRKb6sR4wuzk0Jwwk
YBewt9Tmiwea/TzSL4USiGC+4Dce3o3G0s7js5GKvZYq0f7Tu0EYRUSXxREQ
bXc+khIpgLvKl2Gd3wjYjNqlsr8JEMzpWnIvkcCI/sz2F9EUzgebxrUkE3gi
uCGTx1DwDTIv61X4dppHmelhCmdHqwenFY6eKdNJiqeQ8dVxi8FJApKnIb02
SRT+6svuESqsmvAmJucChddWhYFmKQS4PfVXg+9RKK706jt1SvG3Vpj7+VoK
1eu0/riq8OQ2XWVxHQXVxcHRVQpXVsSF6oopJLL//tansG0az7q9hcLPR1xf
8VMJ7NxwrdXyLYXnwenfBhVeE+nyu+87CgG/+KRJFNYolPPSeyksr3q5XCON
QJfmngsfByiMXWuy3qjwz8NaftljFCaIPzadVRj6Yq2mzxRC7LtrChVmCQ42
fZmksNT+tnmdwoMP20ydpincKz+oM6Rw87cj7w7PUNh13ztRpnAx1zizYI7C
vvLpEc10Aqn+L396IaVwy9/BjaNweHbStFxOIevi1kpbhf8FUT0pYg==
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 1.7000000000000002`},
  PlotRange->{{1, 100}, {1.7014461616150467`, 2.7182754811714283`}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{{3.515497777726169*^9, 3.515497794507268*^9}, {
   3.515498197613374*^9, 3.515498205261358*^9}, 3.515756979260985*^9},
 ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztXAd0FcUa/iHAUwS72PCJBQSkN8tDsNCkSJNeQzOISAtdpBiIdAQB6T2U
0Am9hk5CDwkECKFESAhJSO+5b77Z2b17k527e3M8T945niNJ7u588035+8y1
mbunR+/+7p59e7qXajjYfaBH355DSjUYMJg9citAVMCDiG6VIvxtY3+KHzb2
n60vfuCJ/e/m/Fdx9jOJzvd7n16u8CnVrVuX/atNNWp2p5PxrHXWIzo0oSsN
nDKXJg/oRhP9YyiLoT7i2GfY+zgKWD6SPn+1BNUZu40iM+SYWhxTjGyZ0bSP
vR8wdSEt8O5Lfeeco4RsKU6FZcUcogldutOE+Utp/vju1NXLn2KypLBqHPYi
2bIT6eq+5TSy8otUst0M2nTiEWXa5KOUPC+orFkLdTDRe2hUrwUUkpzDn2TG
BdOd5GxKOt+fKtSaS+FsHTJuzaKaFQfRxWTbx7xRUbKlh9GMWtXJKzSNP5E0
V6ecGNiPKtVZSRGYZkoQDatck6beTJfBamobmk43JlehMj0DKRFTSgyg7mWq
0owwKZSqc+wL7GciBXQuQx8tvacsk3yU0udirVoKGcm4v51G9plMR6Iz+RMM
7+bUKlS6ixhewmnq8F5VmsmGp00hJ54CJ/SgweO9ydtrKDV5601qMHQmbbuf
aI61pdF178pUpsdZSkKbjNs0v3opanY82gI2ix6sb0hlGm+hKCaTWZGbqFHp
RuT7IEWKlS+dbJ7yOYilayXkJf3eeupdpyP9HpREObYUujywLJUdcImSAUu+
SP0+KEeeQSmqvKSHr6Zh80MoJU0vZlKYfdKpFDK6LL1SczoFp+RQxt2VVL/U
pzQ3PF4GtStzViz5j6xD9Qd60YB6n9PoY7GULaeUr5UMEiMdQq61stmyKdL3
MypZ8w+6k8E0YApb4q6OS8w0QBu4LYMiNo+lhZeCaLq2VlKYXkAiVnxMpVv5
kM+SaeTp3ovGbb1NaXKotlZpoV5UvdJQpiM5lHRhMFWs4U3X06Q4Z3JlDEmU
DqFgbpOfQ3FHPajllFBKZU8STrWjdyqMp2tpqlbPo9tcq2drWq2Z/MxI2j1+
NPWpkNuO5cVogsmsQ51PptPNdDaszIfk73OQojKlOBWWdK43lauF7VRe1yjf
ly4kS2H2BUug0+1LU608NizvCGXPc9v72OM0uWcPGj17Pk3q3po8/aKEF4mh
I17u5DljAU0f0oO8j8bCW6j2Pivaj7qVrqytk6S5JpM5qRS2/mfynL6YFrFm
Tfv50v0Mc1j6XVo7qBMNn7OEZnt2oMHr71G6nE3vGYO3zKYR7p2o/56HiveT
w2TPC+aNLPR/v8x/uZESiggG/Gz1v3n/PH7xQKgtfrjpUUUcUezf5UsP6Id+
24oadPqUeHY2MAJNxKe87YobPHv2n3Z/ebt/Od/JU6fu0LChu5ztZExMMjVu
tOyfnXzCd/LkHRo+bHcJgWzZfBU9eJDAP71n0Fvdv/iZa6Mb6rmLjh8Lf0JH
t3DBGVq+7JzyGj9ruQQ/dCiMfhq9N7/wu3cfU4d2ay3D9++/Qb+MP6jCs7Nz
uLomJaW7An9VwPt9v40uXrxvZ6+pgyuL6qbHkQqcNfM4bfQNMgaa8+70u0YT
vQ7lF37rVgx17bzeMnzfvhvkNUGD51k0E/jWLcE0Y/qx1wT8x/7b6dy5P+3s
NVyCz5t7ilavupAHrnhlaW6eHjaLapdvT/5xObK0v5C1ARw5nEtixQCUoAXG
KG+6a8uMo9DL9yglR5pGW2SPjEyk1i1XSdgL5WZ30zRaloRmkQnz64IZRvLh
wyQ7c3WnzGqvxtmcdLE3b75Cs2YcV0mx0AcP3nSV1DgtUl4XFnq54I8zGqsq
2uvXXeK6KadzSANUOklSIZvjmtUXwK7OMSQ4inq4b8xDisUTEp0rptbW1jhK
t8ibkZFNjRospcTEdAdemSIK+Bs6PUZEpA27mkvwRQsD7FtgAb5yxXlavChA
hYNZDZmtwEHFRqDCMeeG9ZfwNXAB/qaA9+65iYKCIu3sVV2Cz/v9FK1aeT6/
cMx94I87jOF6tyNwpAJTUjL4jqelZbrCW1LA+7P1Dgi4Z+et4hw+57cT5Osb
pMIRJcB6W4Uz2wMbpMKDr0SRezff/MLhtZo1Wa6YMOvwtwR8/NgD5Lfjqp29
skvw3btCaezP+/MLj45O5oPHJFyA/1vAe3T3paDLD+zslSQCM23qUdq+LYRU
4B/zT9OSxYHGQANexDgs2FDhSH379N6cXzhm2+KblXz2LsDfFnDviYe5D9PY
K7oE9/e/RUMG7cwvHLrWtPFySkhIcwVeSsDHMWHZtfOanb2CS/DDLKweMXy3
ZfjIEXvoxPHbpXSDb/K1LsQzgbOIHJH5OwI+ZvQ+HjRq7B9K5E3gSAUeYIEq
y7CNgea8qamZfM3j49Ncgb8r4FhDB4Ep7xIcZrlf3635hWdlZXNpj4iItwT3
HLwT1ljNCOcynwKfqrGXcw6H82T6qcLDwmKoTes1dgsj4Dyslh/m6cNqyTGX
LBoRA3hfDKB3r00UGKhzLmV1A7BJT8fskbX08MjiABB5QgJyD6CgeK8/c1K6
RFgtP3eySAv7gGxKk1lBqyhIoVznNaI0CmbpgQg5hLmwwSq1utmY5ob1l+1T
/cCR0/HcQ8cpP5AwmW1pwYzRdO64jnJyHJl57/IDBDfRu6S0bkJeRpAjekCp
VZt2GR25TVqQV8klxWoZOQQ69Fq0Sg63OopZ2dzkMv3EYJlSqnBYZEQAmnWw
Bv9AwBF3IRzT2Eu7BEfW+U3TFXwULsDLCjicmYO8ve8S/Pr1aGrdajWlp2e5
Ai8n4LCzDmb9PR1c740EjlQgIvx2bXy4Zc4DNOBFIerPiPjyAj5owA4e92m8
7zqHt/12DRZahaNGAlXRLLI1+IcCDk908IAuc37HJTi8MsJGTVFdgwPW12OL
MgLr8Aq6uXdov5YyM8XKl7IEryjgqJg67PjbLsFv3ozh1Q4t+jGBN2uygpnv
1EoC/qv3EcfM9t8SeVOLaCpQjfIjHyTmBRrwCnxlAUeG7FCYessleEZGFpc3
ntNbh1cRcPhuKICmpCUtwasKOPKrpUt0mcabLsGjWVbXvNkKunfvsSV4/a8W
Y6TVBBx1H88hulj/DT28cB543c8WiDsaStyGJIdH6xaxNQQWlg3Dhn3lb183
x9bUtvscrwJoGvqaRMwAwu8augF79NmC1DwvUEJaS2Dv3InjQert23EK9lVz
7EcCi4wAJlarPpSQY2H0vqi7kD4WWO9Jh2n6tGP2DXpFjoVKQTxULLKCzp3W
88KtVewnug3CfHHOwFsoIOthcaKsTSEn7J8KdhQcunXdwGfAW7ykY7cWE0va
OGP/j2Cf8qs/r1loAvaiwm4UEBeRPnPOgfobPBUvMug4JNFvYePHzmZTWzDF
xqbwBEeTghccmXLFvIWNH9tD689EvwgYWrVYxau4vN/nxR45j2jtMbW8jdG0
oqIS+TTqCHp4SsRGONPhLZ7T0ZvHtJL3RsywUXCUdXXqjJFoBbVn5WqlYj8X
WBQle/bYaPewxc2xXwgsXKvHd1vs2GLm2C8Fdv68047YZ8yxX8mwRc2x9QQW
BV8H7NPm2PoCOzc39ik5FsIIE9tAYH+fc5JjtRK//mxX9REqCG9U4JzZJ6l7
V1+7YypiTtpIYNf6XORHVCz5UbCFzbFfCyy0Ew7x5AkRhBQyxzYWWBSIoYlb
NgcrWDdzbBOh5vfvJ1AX5iOwUcz38EYFzOFNBRxJ/ADmixF+sX2ymSObCSQM
IE6cOnVYRzeuP8oxRzYX00WJD64JiWW2HCauH4lrZ4V4hATV85pwEOlcpjmy
pWp0WWSHAjHWiZmbDANRUiF400rAYCfUcTJXkmZO+K1AItiA9KLSxWLiVHNk
G/EOIY7Pmot6Xko2h7fTESN7GTzID8RJ5sj24h3ilnUsmAQxMhB8TjCHdxTE
dxkx3GGvHpsgzvHmyE7iHWa4b+91rni4GvDoUTLFyeGI1lHc7iLeIbFHjQDp
x7q1l7hYxsjh4hYGddMtGJIu5A9nTt99ZI50F8ij/uEchcFcu/rwoYFQqRC8
cddt8I7tV7kD+nnMfh6PRpqz9hLvEIIi5keWh3QNZu6+ObyPeAdTjGQFid/M
Gce49Ygwh3uId48fp/JSLSoqXr8covBbsXTXHP69eJecnKHYWOb2R4/ai2Wj
2+bwH/inIjzP28Z0oiPLr7Hqx46GQ0zDlB5cqfeCIO+XNYxiB3BgqD+KMUAt
cCkFeVP7tj5c5JiPuiHGoHhWC+GtGorqvwRR2IBfvX8zSLyDq8BnHBZMZyYt
PDyWlEKN83Kv9CsE6YbBmko7REwbtwNxCgolhW3BEVNmZnaIwiyJed1Ex0Y3
8Z3NdaggRQ0FvgJWFPIONWf6esWRNFecq5IaXWnnpJxzhKBAwQIRPCQSd6JQ
9WIW5ZLYTitRbqHcl8CdTW2U2A+IEcR71Mi9PNafwN4hcsCUz+vI5XFwYbHb
+lvVzpjHCOb4+FTatjWYH4lDjZGbXmImGWY4UK6LajfjRDewO/BTqLQh5sN5
61WmzujmtHk3v4huYEBgj7p12cC7QfSG0WB9Tsi7QXqJJG+S6AZGFFUFXEaB
pGBSOO+HMzgq70bcE6LJohuUqBHgw5HhEBH2ec/uUG72Dpt3M41/eori4lI5
DCko3BK8HdYK1Tm2PvvNe5opeoI4IE/6bdYJbvUwt0kTD/PIhKWFe8x7miN6
wiVRbPuYn/bxqcFPY88CztxDku5n4L3ULvBmLv/5NPdeV65E0rKlZ7kE4U4G
rtPANuCcHu+3mo9qoegO5QHmdXmGglMHqAK0HXsJMUDgttG8u6WiOyQL0CJM
7HuPrXx0GOXiRYHcXSOTXmfe3Ur+qSgf3bmzEXyyMHpfN1zKPT6O5xGbIXNl
grpK3qO4OUQ+okeINVBAoxcUd7Ab6B1+Fd4m+mHSMvMefUWPOHw7feouH+Pw
obu5Y4eRhmHBRR3sL5PHheY9bhU9QgmgplBFLB5mDdWEF1yx/Bw3W8wdzDXv
URGqYnzWiA9xRo40BbU/zBrSjNgL24V30GLWdqZ5z/v4p2e4+T5+/DYvIENZ
sTvY9O+Yd5786xEeKcEOoD49xUDE1f7w5pAYLXY9JCSK2xgcPOEwBIYSHhfr
AduyedMVrpdMoLzMR6sYsuLcieIUb9PGIN4L8i6kjNgzKBGeobSJEUNfYUvH
mHcfKLqHZcedXoRISJ+xtKhr47YY1B2Favg6GA9MENoyXN49dBCKo/jCZ3n0
BQuGe894Pn7cAR4FYScxC7hPRKSQOrQJDY1GuDlQzqBe8lLuijzH7SZOrnYz
26lnwMojLUFig6gTsg77CguB8jnWycOcJlzQQGUwfwg5BBqKiI1FBl7vy8U8
S8TaYTuAwy1O1MhgPdzNaZRo/AUuR3BvMEiQF0g+rC/sHHa8Qb0lnArV+6lT
/LkM791znedHkGuY0vbmdLGCDnkQ6KCeqCPgPaIK+DIYhK++WMRzDjgjeGC8
R7sTTH2wVRDOFuZ0KfzTS9wzIXmAwGHUGD3iUVzdgReGxmAxESODEkPBCvhu
uMxjOuwzYgj008CcVvn0MncIMFCgRdIItzqbKSnu7cDiI47DTCEuCE1gElGB
hzhhE+AF4L9AzYKD2gY2QaO02fRVq1e4GUNyiso6VtnP7yr3e8ijYG7BjxXG
tKEQKHngShD4Ydyh2jhVxNhh6qB/1SxOXWlQgssudhluCfOHgUNYjMQZKRE4
MXecEn3L1gLXPTEGBA1YfqQrUC8claHKhRIoU+wyFoehHOu9pi0F7AH8N7QR
XaPih6FAfSAFquDBr8PXwUbAj6AYgy3BNwmghrB4WFbUiJneODmHgmHjB/L4
VF64AVg+JPIwCNhhOHLsDMImGJFBA/34cLAiBbniFeM7hWoSwg3oA9rCsOAy
DXYJjhN6AUPMkm4nZV5xZVH/NQxF+yFk6AJdwehitjDAEBhcLIDXh2Fu22YN
37B6Xy7iRgjlS5zJOSm1OnAqdd5CXJdgD1FvYxtrFS603pX82fh/WmCUSEqZ
rB0CSb70b5XJ+Q0o2Tfj4w1TYqP+Ta86yb5LnqG76mTUofQek+wL18ZnLoar
b+3CkuRbyi6xWLiZZPwlX6ssRrKNGEQ70XbSTtw7/f9rpzoqrYHqqJo/AYPL
ZzuEBwcstFOv3T4x7dTN0BpY2Iy/axKLFwXYD3sd2v1930inAv8FTyghTw==
\
\>"]]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"L", "[", "n_", "]"}], ":=", 
   RowBox[{
    RowBox[{"F", "[", "Pi", "]"}], "-", 
    RowBox[{"F", "[", 
     RowBox[{"g", "[", "n", "]"}], "]"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.515756858396167*^9, 3.515756882889103*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", 
  RowBox[{"Series", "[", 
   RowBox[{
    RowBox[{"L", "[", "n", "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "Infinity", ",", "3"}], "}"}]}], "]"}], "]"}]], "Input",\

 CellChangeTimes->{{3.515756892491153*^9, 3.515756897323286*^9}, {
  3.515756953495856*^9, 3.515756954275492*^9}, {3.515757013619231*^9, 
  3.515757018301708*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   RowBox[{
    FractionBox["1", "12"], " ", "\[ImaginaryI]", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"3", " ", 
       SuperscriptBox["\[ExponentialE]", "2"]}], "-", 
      RowBox[{"6", " ", "\[ExponentialE]", " ", "\[Pi]"}], "+", 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         SuperscriptBox["\[Pi]", "2"], "-", 
         RowBox[{"6", " ", 
          RowBox[{"PolyLog", "[", 
           RowBox[{"2", ",", 
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{"\[ImaginaryI]", " ", "\[ExponentialE]"}]]}], "]"}]}]}], 
        ")"}]}]}], ")"}]}], "-", 
   FractionBox[
    RowBox[{"\[Pi]", " ", 
     RowBox[{"Log", "[", 
      FractionBox["1", "n"], "]"}], " ", 
     RowBox[{"Log", "[", 
      RowBox[{"2", "-", 
       RowBox[{"2", " ", 
        RowBox[{"Cos", "[", "\[ExponentialE]", "]"}]}]}], "]"}]}], 
    RowBox[{"2", " ", "n"}]], "+", 
   FractionBox[
    RowBox[{
     SuperscriptBox["\[Pi]", "2"], " ", 
     SuperscriptBox[
      RowBox[{"Log", "[", 
       FractionBox["1", "n"], "]"}], "2"], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "\[ExponentialE]"}], " ", 
        RowBox[{"Cot", "[", 
         FractionBox["\[ExponentialE]", "2"], "]"}]}], "+", 
       RowBox[{"Log", "[", 
        RowBox[{"2", "-", 
         RowBox[{"2", " ", 
          RowBox[{"Cos", "[", "\[ExponentialE]", "]"}]}]}], "]"}]}], ")"}]}], 
    
    RowBox[{"4", " ", "\[ExponentialE]", " ", 
     SuperscriptBox["n", "2"]}]], "+", 
   FractionBox[
    RowBox[{
     SuperscriptBox["\[Pi]", "3"], " ", 
     SuperscriptBox[
      RowBox[{"Csc", "[", 
       FractionBox["\[ExponentialE]", "2"], "]"}], "2"], " ", 
     SuperscriptBox[
      RowBox[{"Log", "[", 
       FractionBox["1", "n"], "]"}], "3"], " ", 
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["\[ExponentialE]", "2"], "+", 
       RowBox[{"4", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"Cos", "[", "\[ExponentialE]", "]"}]}], ")"}], " ", 
        RowBox[{"Log", "[", 
         RowBox[{"2", "-", 
          RowBox[{"2", " ", 
           RowBox[{"Cos", "[", "\[ExponentialE]", "]"}]}]}], "]"}]}], "+", 
       RowBox[{"3", " ", "\[ExponentialE]", " ", 
        RowBox[{"Sin", "[", "\[ExponentialE]", "]"}]}]}], ")"}]}], 
    RowBox[{"24", " ", 
     SuperscriptBox["\[ExponentialE]", "2"], " ", 
     SuperscriptBox["n", "3"]}]], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", 
      FractionBox["1", "n"], "]"}], "4"],
    SeriesData[$CellContext`n, 
     DirectedInfinity[1], {}, 0, 4, 1],
    Editable->False]}],
  SeriesData[$CellContext`n, 
   DirectedInfinity[1], {Complex[0, 
      Rational[1, 12]] (3 E^2 - 6 E Pi + 
     2 (Pi^2 - 6 PolyLog[2, E^(Complex[0, 1] E)])), Rational[-1, 2] Pi 
    Log[$CellContext`n^(-1)] Log[2 - 2 Cos[E]], Rational[1, 4] E^(-1) Pi^2 
    Log[$CellContext`n^(-1)]^2 (-E Cot[Rational[1, 2] E] + Log[2 - 2 Cos[E]]),
     Rational[1, 24] E^(-2) Pi^3 Csc[Rational[1, 2] E]^2 
    Log[$CellContext`n^(-1)]^3 (E^2 + 4 (-1 + Cos[E]) Log[2 - 2 Cos[E]] + 
     3 E Sin[E])}, 0, 4, 1],
  Editable->False]], "Output",
 CellChangeTimes->{{3.515756954794751*^9, 3.515757024651775*^9}}]
}, Open  ]]
},
WindowSize->{978, 616},
WindowMargins->{{76, Automatic}, {23, Automatic}},
FrontEndVersion->"8.0 for Linux x86 (32-bit) (February 23, 2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 684, 14, 30, "Input"],
Cell[1266, 38, 838, 24, 50, "Output"]
}, Open  ]],
Cell[2119, 65, 1085, 29, 59, "Input"],
Cell[CellGroupData[{
Cell[3229, 98, 499, 14, 30, "Input"],
Cell[3731, 114, 1152, 33, 63, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4920, 152, 536, 15, 30, "Input"],
Cell[5459, 169, 217, 4, 34, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5713, 178, 356, 7, 30, "Input"],
Cell[6072, 187, 957, 29, 50, "Output"]
}, Open  ]],
Cell[7044, 219, 398, 12, 30, "Input"],
Cell[CellGroupData[{
Cell[7467, 235, 367, 9, 30, "Input"],
Cell[7837, 246, 880, 26, 51, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8754, 277, 615, 18, 58, "Input"],
Cell[9372, 297, 2083, 61, 89, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11492, 363, 708, 20, 62, "Input"],
Cell[12203, 385, 9889, 168, 237, 3348, 59, "CachedBoxData", "BoxData", \
"Output"]
}, Open  ]],
Cell[22107, 556, 280, 8, 30, "Input"],
Cell[CellGroupData[{
Cell[22412, 568, 386, 10, 30, "Input"],
Cell[22801, 580, 3283, 91, 143, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
